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Abstract 

In this paper, we perform a global analysis of i? — > PP and PV decays with 
the QCD factorization approach. It is encouraging to observe that the predictions 
of QCD factorization are in good agreement with experiment. The best fit 7 is 
around 79°. The penguin-to-tree ratio \P-k-k/Tt^-k\ of vr+vr^ decays is preferred to be 
larger than 0.3. We also show the confidence levels for some interesting channels: 

— > vr'^vr'^, K~^K~ and B~^ — > ojtt'^ , ioK~^. For B — > ttK* decays, they are expected 
to have smaller branching ratios with more precise measurements. 
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1 INTRODUCTION 



The charmless two-body B decays play a crucial role in determining the flavor parameters, 
especially the Cabibbo-Kobayashi-Maskawa (CKM) angles 7 and a. With the successful 
running of B factories, many charmless decay channels have been measured with great 
precision. However, since hadronic B decays involve three separate scales, mvi/, ""^b, and 
J^QCDi where perturbative and nonperturbative effects are entangled, it is highly nontrivial 
to relate flavor parameters to experimental observables. 

Recently, theorists have made much progress in nonleptonic B decays: three novel 
methods, QCD factorization (QCDF)pP, the perturbative QCD approach (pQCD)|2] and 
the charming penguin methodjS], have been proposed. These methods have very different 
understandings of B decays: For both the QCDF and pQCD approaches, the factorization 
theorem is proved for nonleptonic B decays in the leading power expansion, i.e., short- 
distance physics related to the scales Mw and rrih can be separated from long-distance 
physics related to the hadronization scale Aqcd, and the long distance part can be param- 
eterized into some universal nonperturbative parameters. In this sense, they are similar. 
But the pQCD approach implements the Sudakov form factor to suppress the end-point 
contributions and proves the factorization theorem in which the form factors are pertur- 
batively calculable. Notice that Sudakov form factor itself is a perturbative quantity; it is 
rather radical and controversial to prove the factorization using the Sudakov form factor, 
while in QCDF the form factors are believed to be nonperturbative parameters. Therefore 
these two methods have completely different power behaviors for B decays. Their pre- 
dictions of B decays are also quite different. For instance, pQCD generally predicts large 
strong phases and direct CP violations, while QCDF favors small direct CP violations 
in general because of the as-suppressed strong phases. The charming penguin process, 
i.e., (bq) (cg)(cs) {<l'(l){(l's), might be potentially important for penguin- dominant 
decays because it is doubly enhanced by CKM factors and Wilson coefficients. The char- 
acteristic of the charming penguin method is that the soft-dominance charming penguin 
plays an indispensable role for penguin- dominant decays. While in QCDF charm penguin 
contributions are hard dominance and therefore perturbatively calculable according to 
naive power counting rules. 

Now BaBar and Belle have accumulated copious data, and will record much more 
data, on nonleptonic B decays. Thus it should be highly interesting to compare the pre- 
dictions of these methods with precise experimental measurements. We gave the QCDF 
predictions on B — PP and PV decays in recent works IHEj. With the experimental data 
at that time, our results prefer a somewhat larger angle 7. For PV decays, the QCDF 
predictions are only marginally consistent with the experimental observation for some 
decay channels. Notice that the QCDF predictions contain large numerical uncertainties 
due to the CKM matrix elements, form factors, and annihilation parameters, and further- 
more, the uncertainties of various decay channels are strongly correlated to each other; 
we are stimulated to do a global analysis in this paper to check the consistency between 
the predictions of QCDF and the updated experimental results. Beneke et al. have 
done a global analysis including tttt, ttK modes with the QCDF approach and have shown 
a satisfactory agreement between the QCDF predictions and experiments, while in this 
work we shall consider not only B PP decays, but also B PV channels. Thereby, 



as we will see later, it leads to some new interesting results. 

One of the most impressive predictions of the QCDF approach is that direct CP 
violation of charmless B decays should be small because the strong interaction phase 
arises solely from radiative corrections. Up to now it has been very consistent with the 
measurements of BaBar and Belle. However, power corrections which may also contribute 
to strong phases are numerically comparable with radiative corrections. Notice that the 
power corrections are difficult to estimate because they generally break factorization. It 
means that the predictions of QCDF on direct CP violations are probably qualitative. 
Therefore in this paper we will not consider experimental results on direct CP violations. 

Our global fit shows that QCDF has an excellent performance on B ^ PP (two 
light pseudoscalars) decays except for the channel rjK^ . But we do not worry 

about it because of the hard-to-estimate contributions from the digluon mechanism and 
the potential large power corrections in this channel. The CKM angle 7 is preferred to 
be around 79° which is slightly larger but still consistent with the standard CKM global 
analysisjH]. We also discuss the preferred range of the penguin-to-tree ratio \P-k-k/T.„t,\ 
^ which is crucial for the extraction of angle a. For B PV (one pseudoscalar, one 
vector) decays, QCDF has also a good performance where the annihilation topology plays 
an important role especially for penguin-dominated decays. But for B nK* channels, 
the QCDF results seem smaller compared with the experimental measurements. How- 
ever, presently there are large experimental errors on these channels, so it would be very 
interesting for BaBar and Belle to update their measurements with higher precision on 
these decay modes. Based on the global fit, we also give the confidence levels for some 
interesting decay channels: 5° vr^vr^, K~^K~ , and B~^ ujtt^, ujK^ . 

This paper is organized as follows: in Sec. H, we will first recapitulate the mainpoint 
of QCD factorization for charmless two-body B decays. In Sec. Ill, the relevant input 
parameters are discussed. Then the numerical results of the global fit and brief remarks 
are presented in Sec. IV. Section V is devoted to the conclusions. 



2 QCD FACTORIZATION FOR CHARMLESS B 
DECAYS 

As we know, charmless B decays contain three distinct scales: M\y ^ S> ^qcd- To go 
beyond the naive model estimation, it is important to show that the physics of different 
scales can be separated from each other. This process is generally called "factorization" . 

It is well known that, with the help of the operator product expansion and renormal- 
ization group equation, the effective Lagrangian can be obtained, in which short-distance 

^The defination of the penguin-to-tree ratio \P^t^/Tt^^\ for B^^tt+tt^ decay is |7] 

P.. _ -1 {a% + rial) + (g^o + r^gj) + r^jbg + 264 - ^^f ^ - 6f ^)] 

T.x Rh (fli + < + rial) + «o + rial) + rA[bi + 63 + 2&4 - \{h^^ - bf^)\ ' 

where Rb = = VP^TW, and rA^-M^- 



effects involving large virtual momenta of the loop corrections from the scale M^y down 
to fi = 0{mh) are cleanly integrated into the Wilson coefficients. Then the amplitude for 
the decay B M1M2 can be expressed as 

A{B^M,M,) = E Ka{fi){M,M2mfi)\B), (1) 

V i q=u,c 

where A, is a CKM factor, Cj(/i) is the Wilson coefficient which is perturbatively calculable 
from first principles, and {MiM2\Qi{^)\B) is a hadronic matrix element which contains 
physics from the scale /i = 0{mij) down to Aqcd- In a sense, this process may be called 
"first step factorization" . But it is still highly nontrivial to estimate the hadronic matrix 
elements reliably because the perturbative and nonperturbative effects related to nib and 
Aqcd are strongly entangled. 

Three years ago, Beneke et al. put forward the QCDF approach in the heavy quark 
limit for B ^ titt \T\. They show that, neglecting power corrections in 1/m;,, the hadronic 
matrix elements can be factorized into hard radiative corrections and a nonperturbative 
part parameterized by the form factors and meson light cone distribution amplitudes. In 
the following we will outline their reasoning. 

2.1 QCDF in the heavy quark hmit 

First, we need to have some knowledge about the end-point behavior of the light cone 
distribution amplitudes (LCDAs) of the mesons. At the scale of rrib, the LCDAs of the final 
light mesons — for example (f){x) of vr or mesons — should be similar to the asymptotic 
form. Therefore it is reasonable to assume that the end-point of the LCDAs of the light 
mesons is suppressed by A/m^. For B mesons, the spectator quark is assumed to be soft 
and have no hard tail, i.e., (f){^) ~ mt/A, for ^ < A/mf, and </)(^) = 0, for ^ > A/rrib- 
With the above assumptions, the form factor is argued to be nonperturbative dominant; 
thereafter, naive power counting rules are constructed and the leading power radiative 
contributions in l/nib can be identified (see Fig. 1). 

Notice that, in Fig. 1, the emission meson from the decay vertex carries large energy 
and momentum (about 771^/2) and therefore can be described by leading twist-2 LCDA in 
the leading power approximation. For factorization to be held, these radiative contribu- 
tions should be hard dominant. For vertex corrections (Figs. 1(a)- 1(d)), every individual 
diagram contains infrared divergence, but these infrared divergences are canceled after 
summation. This cancellation is not accidental. Intuitively, the qq pair of the energenic 
emission meson can be viewed as a small color dipole. Since soft gluons can not taste the 
difference between a small color dipole and a color singlet, the emission meson decouples 
with the soft gluon interaction. This argument is well known as "color transparency" 
[TH] . Technically, not only soft divergence but also coUinear divergence is canceled. For 
penguin corrections (Figs. 1(e) and 1(f)) and hard spectator scattering (Figs. 1(g) and 
1(h)), since the end point of the twist-2 LCDA of the light meson is A/ruf, suppressed, it 
is not difficult to show hard dominance. So factorization does hold in the heavy quark 
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Figure 1: Order of corrections to the hard scattering kernels. The upward quark 
hues represent the emission meson from the b quark decay vertex. These diagrams are 
commonly called vertex corrections, penguin corrections and hard spectator scattering 
diagrams for (a)-(d), (e) and (f), and (g) and (h), respectively. 

limit, and the corresponding formula can be explicitly expressed as 

1 1 
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+ J didxdyTl\i, X, y)^B{O^MAx)^M.Ay) 



= (M1M2IJ1® J2|5)-[l + 5]r„< + 0(AQCDM)]. (2) 

In the above formula, ^b{C) ^i-nd (^Mi{x){i = 1,2) are the leading twist wave functions of 
B and the light mesons, respectively, and T/'^^ denote hard scattering kernels which are 
perturbatively calculable. The readers may refer to Ref.^T] for more details. 

One of the most interesting results of the QCDF approach is that, in the heavy quark 
limit, strong phases are short dominant and arise solely from vertex and penguin correc- 
tions which are at the order of a^. It means that, for charmless hadronic decays, direct CP 
violations are generally small because strong phases are as suppressed compared to the 
leading "naive factorization" contributions. But in principle power corrections may also 
contribute to strong phases, and numerically AqcD/fnb is comparable to Furthermore, 
there is no known systematic way to estimate power suppressed contributions (note that 
soft collinear effective theory ^2] may be a potential tool), so QCDF could only predict 
strong phases qualitatively. 



2.2 Chirally enhanced power corrections 

The above discussions are based on the heavy quark hmit; i.e., power corrections in 1/mb 
are assumed to be neghgible. Then the question is, for phenomenological apphcations, 
whether it is a good approximation. There are various sources which may contribute to 
power corrections in l/rrib] examples are higher twist distribution amphtudes, transverse 
momenta of quarks in the hght meson, annihilation diagrams, etc. At first sight, power 
corrections seem really small because they are suppressed by AgcD/fnb — 1/15. However, 
this is not true. For instance, the contributions of operator Qq to decay amplitudes would 
formally vanish in the strict heavy quark limit. But it is numerically very important 
in penguin-dominated B rare decays, such as the interesting B irK decays. This is 
because Qq is always multiplied by a formally power-suppressed but chirally enhanced 
factor = ^^j-^^:^^^) ~ C^(l)) where mi and m2 are current quark masses. Another 
example is annihilation topology (Fig. 2), the importance of which was noticed first in 
the pQCD method 0. Therefore phenomenological applicability of QCD factorization in 
B rare decays requires at least a consistent inclusion of chirally enhanced corrections and 
annihilation contributions. 




Figure 2: Order of ag corrections to the weak annihilations. 



Chirally enhanced corrections arise from twist-3 light cone distribution amplitudes; 
thus, the final light mesons should be described by leading twist and twist-3 distribution 
amplitudes. Then we need to redemonstrate that the leading power radiative corrections 
(Fig. 1) are still dominated by hard gluon exchange. Unfortunately it is not true for hard 
spectator scattering which contains logarithmic divergence in the end-point region. A 
similar divergence also appears in the annihilation contributions. It means that, strictly 
speaking, factorization does not hold for chirally enhanced corrections and annihilation 
topology. The readers may refer to Refs. [71ll3j for more technical details. Phenomenolog- 
ically, Beneke et al. [3| introduced a model parametrization for the end-point divergence: 

J X Ah 
" 

where Xa denotes the annihilation contribution and Xh denotes hard spectator scattering. 
We will follow their approach in this paper. 

For the rest of the power corrections, they are argued to be generally small ^3] based 
on a model estimation with renormalon calculus. 



With the above discussions, the decay amphtudes can be written as 



10 

A{B ^ M1M2) = ^ E J2M(^'iMiM2mB)f + IbImJuA). (4) 

V p=u,c i=l 

where {MiM2\0i\B) f is the factorized hadronic matrix element which has the same def- 
inition as that in the naive factorization approach. For the exphcit expressions of QCD 
coefficients and annihilation parameters bi, the readers may refer to Refs. |31EIIZ| ^- 



3 INPUT PARAMETERS 

The decay amplitude for B M1M2 depends on various parameters, such as the CKM 
matrix elements, decay constants, form factors, renormalization scale /i, LCDAs, and so 
on. Notice that although the predictions of QCDF are formally scale independent at 
one-loop order, numerically there is still a small residual dependence. In the global fit, 
the scale /x is varied from mh/2 to 2mb. For the rest of the parameters, we will specify 
them in the following. 



3.1 CKM matrix elements 

The CKM matrix in the Wolfenstein parametrization is read as 



Vc 



CKM 



( I-AV2 A A\'^{p-i7]) 

-A 1 - AV2 I + O(A^). (5) 

\ AX'^il- p-iT]) -A\^ 1 



It contains four parameters A, A, p and 77, in which the first two are well determined [T^ : 

A = IK.I = 0.2200 ± 0.0025, AX^ = \Vcb\ = (40.4 ± 1.3,*,* ± 0.9*^,^0) x 10^1 

As to p and i], they are kept free except for the constraint \Vub\ = (3.49 ± 0.24stat ± 
0-55theo) X 10-3 [16J and sin 2/? = 0.731 ± 0.055 (world averaged) |8j. 



3.2 Form factors and decay constants 

The form factors and decay constants are nonperturbative parameters. The form factors 
can be extracted from the semileptonic decays and/or estimated with some well-defined 
theories, such as lattice calculations, QCD sum rules, etc. But the related errors are still 
sizable. The decay constants can be extracted from the leptonic or electromagnetic decay 
width with high precision. In the fit, we choose the corresponding numerical values as 

follows [niiiHiiiniiini: 

= 131 MeV, fK = 160 MeV, Jk* = 214 MeV, 

^In fact, there is minor difference for the hard spectator scattering term /^^ between ref.jjj and ref.|7j. 
This is due to a subtle point relating to some divergent boundary terms in a process of integration by 
part In this paper, we adopt the expression of /^^ from ref.[7]. 



fp = 210 MeV, = 195 MeV, ^ = 233 MeV, 

fg = 140 MeV, /, = 176 MeV, = 39.3°, 
/s = (180 ± 40) MeV, R^k = 0.9 ±0.1, 

F^J{0) = 0.28±0.05, ^^^'(0) = 0.47±0.07, 

= 0.37±0.06, ^^"(0) = ^^''(0), 



where R^k = f^F^^/fRF^^- In the above, we assume ideal mixine; between u and (h, 
i.e., uj = {uu + dd)/\/2 and (f) = ss. As for rj — rj' mixing, we follow the convention in the 
quark-flavor basis fIT\ |22] and assume that the charm quark content in 7]^'^ is negligible, 

{0\q7til5q\v'''Kp)) = ifiioPt^ = s). 



2 



" ^^^^^^ " ^ 

pB^/cos^ sin^N o^/ ^^/sin^ cos^N 

where the four octet-singlet parameters can be related to three quark-flavor parameters: 



fs = ^l/3f^ + 2/3/2, /o = y'2/3/2 + 1/3/2, 

6*8 = - arctan {V2fs/ fq), Oq = (p - arctan iV2fJ /,). (7) 

3.3 LCDAs of the mesons 

The LCDAs of the mesons are basic input parameters in the QCDF approach. The LCDAs 
of a light pseudoscalar meson are deflned as 



{P{p')\UyMx)\o) : 

4 JO 



1^ (ly^ ^i(up' -y+up' -x) 



X S/750(m) - /ip75 - (y^vP^^z""- 
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(5q 

.2 



where z = y — x, (j){u) {(f)p^^{u)) is leading twist (twist-3) LCDA, and hp = rrip/ (mi(/i) + 
7712 (/i)) (here mi and m2 are current masses of the valence quarks of the pseudoscalar 
meson). Because the current masses of light quarks are difficult to flx, we would like to 
take 



-Cu 
J_q 



which is numerically a good approximation. For the related quark masses, we shall follow 
Ref. W\: 

m,(2 GeV) = (110 ± 25) MeV, mdmb) = 1.3 GeV, mfe(mb) = 4.2 GeV. 



For vector mesons, only longitudinal polarization is involved in B PV decays. 
Furthermore, the contributions of twist-3 LCDAs of vector mesons are doubly suppressed 
by as and A/rrib', therefore, they can be safely disregarded. Then the leading twist LCDA 
of a longitudinal vector meson is defined as 12^ 

mp')\UyMx)\o) = ^ £ dne^^-^'-y^^^'-%\{u)^UE. (9) 

We shall use the asymptotic forms of the LCDAs for the following discussions: 

= (p\\{u) = 6uu, 0p(m) = 1, (pa{u) = 6uu. (10) 

Strictly speaking, the asymptotic forms are only valid for /i — > oo. We notice that in 
Ref. ilj, Beneke et al. employ an expansion in Gegenbauer polynomials for leading twist 
TT, K LCDAs. However, since there are many light mesons involved in our global fit, if we 
consider a similar expansion for the leading twist LCDAs, many free parameters would be 
introduced. Fortunately, the corrections to the asymptotic form are numerically not so 
important because they only affect part of the vertex and penguin corrections (numerically, 
the readers may refer to Table 3 of Ref. [Tj to see the effects of the Gegenbauer expansion). 
So for simplification, only the asymptotic forms are used in our discussions. 

For the wave function of the B meson, only the moment /q^ (i^^s (0 = "^b/^b 
appears in the factorization formulas. We do not know much about the parameter 
and the estimation of Ref. [7j is quoted : = (350 ± 150) MeV. 

4 GLOBAL ANALYSIS OF CHARMLESS B DE- 
CAYS 

In this work, the global analysis is based on the CKMFitter package ^ developed by 
Hocker et al. [B]. The original package includes B hh {h = ir or K) decay channels, 
and we enlarge it to include B PV and B riTT{K) decay modes. The Rfit scheme 
is implemented for statistical treatment. Simply speaking, the Rfit scheme assumes the 
experimental errors to be pure Gaussians (if the systematic errors are not so large) while 
the theoretical parameters vary freely in a given range. In this spirit, it is similar to the 
95% scan method. One of the main differences is that the overall Xmin is assumed to be 
Gaussian distributed in the 95% scan method, while for the Rfit scheme, the confidence 
level of the overal Xmin computed by means of a Monte Carlo simulation. The readers 
may refer to Ref. [0] for details about the Rfit scheme. As to the QCDF expressions for 
the related decay amplitudes, the readers may refer to Refs. jH E] for B ^ PP decays 
and P for B — i> PV decays. 

Compared with pQCD, QCDF requires more input parameters, such as form factors, 
annihilation parameter Xa, hard spectator parameter Xh, and so on. To make the global 
analysis appear more persuasive and at the same time save computing time, we minimize 
the number of variables by fixing those insensitive parameters. One example is that we 

3http://ckmfitter.in2p3.fr 



shall use the asymptotic forms of the LCDAs for all final light mesons, and not employ 
an expansion in Gegenbauer polynomials. 

Since power corrections violate factorization, the parameters Xa and Xh are intro- 
duced as a model parameterization. But we should be care that, in principle, these 
parameters are channel dependent. Fortunately, assuming "factorized" SU(3) breaking, 
we can see that X^ and X^ are universal separately for B —>■ PP and B —>■ PV decay 
modes. However, there is no way to relate the chiral parameters of the PV channels to 
those of the PP channels. So we have to introduce, besides X^^ and Xf^^, the additional 
parameters X^^ and X^^ for B PV decays. 

In Refs. 13 I2H]; Beneke et al. present a detailed analysis of -B ^ tttt, itK with 
the QCDF approach. They show an impressive agreement between experiments and the 
QCDF predictions: ~ 0-5 for six decay channels. Their best fit results favor 7 around 
90° which seems not so consistent with the standard global fit of the CKM matrix elements 
using information from semileptonic B decays, K-K mixing and B-B mixing. Even for 
7 around 60°, ^ 1 is still good enough to be acceptable. In QCDF, these six decay 
channels are sensitive to several input parameters: the CKM parameters \Vub\ and angle 
7 (or equivalently p and 17), form factors F^"^ and F^^ , annihilation-related parameters 
Xa and Jb/^b, and current quark mass m^. These parameters vary freely only in a given 
range which is either determined by experimental measurements (|Vub|), estimations with 
QCD sum rules and/or lattice calculations (form factors and decay constants), or well- 
educated guesswork (Xa)- So it is really nontrivial for the achieved agreement between 
the QCDF predictions and the experimental measurements. 

In this work, we will extend the global analysis to include 14 B ^ PP and PV 
decay modes (see Table 1). Notice that the hard spectator parameter Xh is numerically 
unimportant for the branching ratios except for a2-related tree-dominated decays and 
even for a2-related decays, it brings at most 20% uncertainties. So compared with the 
global analysis oi B ^ tttt, tcK [7 , our extension would include seven B PV channels 
and newly observed B ^ rjT^ decay, while only three new sensitive parameters — the form 
factor Aq'^ and complex variable X'^ — would be involved. Therefore we can have a more 
stringent test of the QCDF predictions which should give some interesting information. 

Recently BaBar and Belle also gave a strong constraint on direct CP violations for 
many charmless hadronic B decay channels. Their search show that direct CP-violating 
asymmetries are generally small. Within the QCDF framework, strong phases are ei- 
ther Us or A/mft suppressed, which also lead to small direct CP violations in general. 
However, only radiative corrections are perturbatively computable in QCDF, while power 
corrections break the factorization in general. Considering that A/m^ is numerically com- 
parable with as, the QCDF calculations on direct CP violations are probably qualitative. 
So the experimental constraints on direct CP violations are not implemented in this global 
analysis. 

Before doing the global fit, the readers may notice that some decay modes are not 
included in the global analysis although they have been observed. These decay channels 
are listed in Table 2, and we will discuss these channels later. 



Table 1: Experimental data of CP-averaged branching ratios for some charmless B decay 
modes in units of 10~^. The following decay channels are the experimental input of the 
global fit. 

BF(xlO^) CLEO [23] BaBar [26j Belle Average~ 

4.3li:^±0.5 4.7 ±0.6 ±0.2 5.4 ± 1.2 ±0.5 4.77 ±0.54 

5.4lii±1.5 5.5tJ:^±0.6 7.4l|^±0.9 5.78 ±0.95 

17.2111 ±1.2 17.9 ±0.9 ±0.7 22.5 ± 1.9 ± 1.8 18.5 ± 1.0 

11.6^1? t\i I2.8t\i ± 1-0 IS.Otli ± 1-3 12.7 ± 1.2 

18.2tto±l-6 17.5tJ:^±1.3 19.4l|J±1.6 18.1 ± 1.7 

U.Gtlitlt 10.4 ±1.5 ±0.8 8.0l|?±1.6 10.2 ± 1.5 

< 5.7 ' < 5.2 5.3t?;?(< 8.2) < 5.2 
27.6t?:H4.2 28.9 ±5.4 ±4.3 20.811^ 25.4 ± 4.3 
10.4t|^±2.1 < 39 8.0l|jj±0.7 8.6 ± 2.0 
16.0;^;^ ± 2.8 ll-Stle -H 13-1 ± 4.7 
5.5;?] ±0.6 9.2 ±1.0 ±0.8 10.7±1.0;?| 8.9 ±1.0 
5.4l|^±0.7 8.71}:^ ±0.9 lO.OlJ;? 8.6 ±1.3 

< 10 <6.2 <6.2 
<21 5.9li:^±0.9 5.9 ±1.9 
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4.1 Main results of the global fit 

When the decay channels in Table 1 are concerned, the global fit shows that the QCDF 
predictions are well consistent with the experimental measurements: The results in the 
(p, 77) plane are shown in Fig. 3 where Xmin = 4.2 for 14 decay channels. As an illustration, 
in Table 3, we list the best fit values of the global analysis for the related B — > PP,PV 
decay modes with and without chiral-related contributions. Notice that two sets of best 
fit values (with or without chirally enhanced contributions) are obtained with different 
input parameters. It indicates that the newly observed B^ —>■ ujK^ decay can be included 
in the global fit without any difficulty, and that it is hopeful that the decay B^ rjp'^ will 
be observed soon. The corresponding theoretical inputs for the best fit values including 
chirally enhanced corrections are also reasonable: \Vub\ = 3.57 x 10~^, 7 = 79°, F^'" = 



Table 2: Measurements which are not included in the global analysis. 

BF(xlO'^) CLEO 25 BaBar [26j Belle [27j Average 

B+ r]K+ < 6.9 < 6.4 5.2t\i{< 7.7) < 6.4 

B+^7r+K*° < 16 15.5 ±3.4 ±1.8 16.21^;^ ± 2.4 15.8 ± 3.0 

B^^TT-K*+ 161^ ±2 26.0 ±8.3 ±3.5 19.0 ± 4.9 

B+^r]K*+ 26.4^1^ ±3.3 22.11^^ ± 3.3 26.51?;^ ± 3.0 25.4 ± 5.3 

B^^rjK*^ 13.8ti^±1.6 19.8t|^±1.7 16.51^2 ± 1-2 16.4 ± 3.0 

B+ ujK+ < 8 < 4 9.2l|^ ± 1.0 

5+^cj7r+ 11.3;|^±1.5 6.6l?:^±0.7 < 8.2 



which means 



0.24, A^" = 0.31, rris = 85 MeV, /i = 2.5 GeV, fe = 220 MeV, p^^ = 0.5, (f)^^ = 10°, 
= 1, = —30°. As to F^^ , there is no strong constraint and the range [0.24, 0.30] 
is acceptable from the current global analysis. Since the hard spectator contributions are 
nearly negligible except for a2-related tree-dominant decays {B — >■ tt+tt^, 7r~^p^), and even 
for a2-related tree decays, it brings at most 20% uncertainties to the branching ratios 0, 
the global analysis can not give a strong constraint on the hard spectator parameter X^- 
In principle, the chirally enhanced corrections could lead to large strong phases from the 
imaginary part of the annihilation topologies. But since the best fit parameters show a 
small imaginary part: = 10°, = —30°; the global analysis still prefers small 

direct CP violations, which is consistent with the current experimental observations. In 
Ref . [23 , it is argued that chirally enhanced corrections are not indispensable for B ^ Tin 
and 71 K decays. However, we can see from Table 3 that, especially for penguin-dominated 
B — »• PV decays, chirally enhanced contributions play an important role. Note that this 
point is not firmly established: There are significant experimental errors in B^ —>■ K~^p~ , 
ujK^ decays. If these two channels were excluded, we could see from Table 3 that it is 
still acceptable without chiral-related contributions. However, t^^K*^ decay also implies 
large chiral contributions: Without chirally enhanced corrections, it is clear that 

B{B+ 71+ K*^) ~ ^B{B+ 71+ K^) ~ 5 X lO"*^ . (12) 

It is 3 times smaller than the experimentally central value (see Table 2). Thereby further 
measurements with higher precision on the penguin-dominated B —>■ PV decays will 
clarify the role of the chirally enhanced contributions. We will return back to 7cK* channel 
later and explain why we do not include this mode in the global fit. 

It is known that the penguin-to-tree ratio \P„^/Ttj^\ is very useful for extraction of the 
CKM angle a In Ref. [H, the authors show that \P^^/T^^\ = (28.5±5.1 ^5.7)% with 
the QCDF approach using the default values for the chirally enhanced corrections, i.e., 
Xh = Xa = ln(mB/A/i). When considering the uncertainties of the X^ h parameters, the 
theorectical errors would be even larger. So it should be very interesting to obtain the 
preferred \P/T\ ratio ^ from the global analysis. Until now the asymptotic LCDAs were 
used for the global fit because the branching ratios are numerically not so sensitive to the 
corrections to the asymptotic form. But for the penguin-to-tree ratio, the case is different 
and we should consider the Gegenbauer polynomial expansion for the leading twist LCDAs 
of the TT, K mesons. We find that, compared with the estimation \P/T\ = 0.276 ± 0.064 
j2n] (including SU{3) breaking effects), the global fit prefers a suprisingly large value: 
\P^^/Tj^^\ ^ 0.41. The reason may be that penguin annihilation effects increase the 
penguin amplitudes, as discussed in [7] . Considering the relatively large model dependence 

*The decay amplitudes fot i?^ tt+tt" are 



where 6t and dp are strong phases. 



Table 3: The best fit values using the global analysis with and without chiral-related 
contributions for B PP and PV decays. "No chiral" means the best fit value neglecting 
the chirally enhanced hard spectator contributions and the annihilation topology. The 
branching ratios are in units of 10~^. The experimental data are the uncorrelated average 
of measurements of BaBar, Belle, and CLEO (see the data in the last column of Table 1). 



Mode 


50 _^ ^+^- 


B+ ^ tt+ttO 


50 ^ 


B+ ^ ir+TfO 


B+ ^ i^o^+ 


Expt. 


4.77 ±0.54 


5.78 ±0.95 


18.5 ±1.0 


12.7 ± 1.2 


18.1 ± 1.7 


Best fit 


4.82 


5.35 


19.0 


11.4 


20.1 


No chiral 


5.68 


3.25 


18.8 


12.6 


20.2 


Mode 




B+ r]7T+ 




B+ ^ p07r+ 


B+ T]p+ 


Expt. 


10.2 ± 1.5 


< 5.2 


25.4 ±4.3 


8.6 ±2.0 


< 6.2 


Best fit 


8.2 


2.8 


26.7 


8.9 


4.6 


No chiral 


7.3 


1.8 


29.5 


8.5 


3.8 


Mode 


B+ (f)K+ 




5° ^ K+p- 






Expt. 


8.9 ± 1.0 


8.6 ±1.3 


13.1 ±4.7 


5.9 ± 1.9 




Best fit 


8.9 


8.4 


12.1 


6.3 




No chiral 


7.1 


6.7 


5.1 


1.2 





of this ratio within the QCDF framework, the best fit value may be not so meaningful 
to reduce the ambiguity in the determination of sin 2a. However, even assuming that it 
is acceptable for the global fit with x^/^dof < 1 (^do/ denotes the number of degrees), 
the ratio |P^^/T^^| is still larger than 0.3. This result is quite interesting, although 
undoubtedly it needs further tests with larger data samples. 

The confidence levels of the angle 7 and some interesting decay channels are given in 
Fig. 3 and 4, respectively. It is encouraging that the favored angle 7 is around 79° which 
is somewhat larger but still consistent with the standard CKM global fit. B{B^ tt^h^), 
which is crucial for a clean extraction of the angle a, is predicted to be around 1 x 10~^. 
So it is hopeful to be observed in the near future. For the pure annihilation decay 
B^ K~^K~ , although the updated upper limit has been very stringent, the branching 
ratio is predicted to be still several times smaller than that: roughly 10~^. 

It is interesting to have a brief look at the relevant PQCD results. Based on tttt, ttK 
decays, PQCD could extract the central values of the CKM angles |30j: a = 78°, /3 = 26°, 
7 = 76°. These results are (probably coincident) consistent with our best fit results. But 
B{B° 7r°7r°) is about 0.3 x 10"^ in the PQCD method, which is smaller than the QCDF 
prediction and could be tested soon in the near future. PQCD also predicts a large direct 
CP violation for tt+tt^ decay. 




7 [degrees] 



Figure 3: Left plot: confidence level in the (p, ?7) plane for the global fit. The contours 
with shaded area inside indicate the regions of > 90% and > 5% Confidence levels. Right 
plot: the confidence level for angle 7 in units of degree. 






Figure 4: The confidence levels for some selected decay channels. The branching ratios 
are in units of 10~^. The dashed lines denote the experimental upper limits; the gray 
bands denote the experimental measurements with la error. 



4.2 Decay modes listed in Table 2 



Now let us discuss in detail why we do not include the decay modes listed in Table 2 for 
the global analysis. 

For B'^ ujTT^{K~^), the measurements of Belle are not so consistent with those of 
BaBar and CLEO. Assuming A^'^ = Aq'' which should be a good approximation, the 
confidence levels for these two decay channels are shown in Fig. 4. For ujtt^ decay, the 
best fit value is 6.66 x 10^^ which is consistent with the measurements of BaBar and 
Belle. The best fit value for uK^ decay is 6.25 x 10^^ which is consistent with the Belle 
observation but larger than the upper limit given by BaBar. 

As to -B — >■ r]K* decays, the branching ratios depend on the form factor Aq^* . Note 
that there are no other observed PV decay channels relying on this form factor; it is more 
or less trivial to include B rjK* decays, since in some sense Aq^* acts essentially as 
a free parameter in the global fit especially considering the large experimental errors in 
these decay channels. One might argue that the observed B (pK* decay also depends 

But this decay channel requires a new annihilation parameter 
VV decays. Hence it may be better to first have a restricted constraint 
VV decay modes to be observed. Then more precise data on 
B (f)K*, rjK* could overconstrain the form factor Aq^* and give a more stringent test 
of the QCDF approach. 

The decay mode B^ — ^ rjK^ was recently observed by Belle. Within the QCDF 
framework, the corresponding amplitude is proportional to /i^-F'^^(a4+r^a6)+/^-F^^(a4+ 
r^Og). Notice that < 0, and large cancellation occurs which leads to a much smaller 
branching ratio compared with the experimental data. But we do not worry about it due 
to several reasons. 



on the form factor Aq^ 
XY for B 
on with more B 



First, it relates to the special property of rj'^''^ which has anomaly coupling to two 
gluons. Specifically, the digluon fusion mechanism where one gluon comes from 
the b s decay vertex and the other from the spectator quark, is presumed to 
account for the large branching ratios of — > rj'K decays. Although it is arguable 
whether the digluon mechanism is perturbatively calculable or not, the contribution 
should be proportional to the coupling {0\GG\r]^'^) . As we know that 



I 4tt 



GG\ri) 



Mr, 



cot I 



1 

5 



and 



^expt(^ 



1 

4' 



this means that 



A^{B riK) 



A^^P^B^riK) 



A^iB ri'K) 



where A^ denotes the amplitude of the digluon fusion mechanism. So if the digluon 
mechanism were important for B tj'K, it should be also important for B rjK. 



Second, we know that when the leading power terms are abnormally small, the 
next-to- leading power contributions become potentially important. Remembering 
that there is no known systematic way to estimate power corrections, the QCDF 
estimation of this channel is probably correct only at the order of magnitude. 



• Recently, Beneke j2Hl propose a novel possibility: for the annihilation contributions, 
two gluons may radiate from the spectator quark and form a i]^'^ meson. In this 
case, it is the leading power contribution. Furthermore, it breaks the factoriza- 
tion and therefore a new nonperturbative parameter is needed to parameterize its 
contribution. 

From the above discussions, it is clear that theoretically great efforts are needed to quan- 
titatively understand B rj^'^K decays. 

The real trouble is -B ^ ttK* decays. Let us take n^K*^ as an example. In 

QCDF, approximately we have 

A{B+ ^ n+K*') cx ajK'F''^ x const + fBfx'fMV, P) , (13) 
AiB+ cx a4/0F^^ X const + fnUfMV, P) . (14) 

Assuming F^^ / F^^ ^ /i^//^, then 

A{B+ n+K*^) fK*F^^ 

So B{B^ —>■ n^K*^) should be smaller than or at most comparable with B{B^ (j)K^). 
Unfortunately, the updated experimental measurements do not support it: 

B{B+ 7T+K*°) = (15.8 ± 3.0)xl0-^ , B{B+ (j)K+) = (8.8 ± 1.0)xlO-^ 

General speaking, we do not anticipate any novel mechanism for this channel because it 
would have a similar influence on B ^ ttK decays. So presently there is nothing we can do 
from the QCDF side. With the global fit, the confidence level for B ttK* is displayed in 
Fig. 4, from which we can see that B{B~^ 7t~^K*^) is comparable with B{B'^ (pK^) 
as expected. The preferred branching ratios are somewhat smaller than the experimental 
measurements for both ir'^K*^ and 7t~K*~^ modes. Fortunately the current experimental 
errors are quite large. We anticipate that further precise measurements would prefer 
smaller branching ratios for these two decay channels. 

It is interesting to notice that, although there are essential difference between the 
PQCD and QCDF methods, PQCD also predicts that B{B+ ti+K*^) is comparable 
with B{B+ (j)K+) which is about 10 x 10"^. 

5 SUMMARY 

QCD factorization is a promising method for charmless two-body B decays, which are 
crucial for the determination of the unitarity triangle. With the successful running of B 
factories, many B — >■ PP and PV decay modes have been observed. We can do a global 
analysis to check whether the predictions of the QCD factorization are consistent with the 
measurements. Since there are many parameters involved in the global analysis, we try 
to minimize the number of free QCD parameters to make the global analysis appear more 
persuasive. Hence, the asymptotic forms are used for the light cone distribution ampli- 
tudes of the light pseudoscalar and vector mesons. For chirally enhanced parameters, it 



is a good approximation to take = r\ = — '^nr^ / mb{ms + rUy). Assuming "factor- 
ized" SU(3) breaking, the chirally enhanced parameters X^.h are separately universal for 
B — > PP and PV decays. However, and X^jj are independent parameters because 

there is no approximate symmetry to relate B — > PP and B — > PV decays. 

With the above set of parameters, we enlarged the CKMFitter package to include 
more charmless decay channels and did a global analysis. It is shown that the QCDF 
predictions are basically in good agreement with the experiments. It is encouraging to 
see that the favored angle 7 is roughly consistent with the standard CKM global fit. It is 
quite interesting to see that the chirally enhanced corrections may play an important role 
in penguin-dominated B — > PV decays. The penguin-to-tree ratio iPTrTr/T^Trl is important 
for the extraction of the CKM angle a. The global analysis favors this ratio to be larger 
than 0.3 with per number of degree smaller than 1. 

Notice that the observed decays r]K^, r]K*, ttK*, ujtt^ , ujK^ are not included in the 
analysis. In this paper we discussed these decay channels in detail. Among these channels, 
only the decays B nK* are somewhat troublesome. In QCDF, B{B~^ tt+X*") 
should be smaller than or at most comparable with B{B'^ — > (t)K'^), which is not so 
consistent with experimental observations. In fact, the global analysis prefers somewhat 
smaller branching ratios for both ti^K*^ and ti~K*^ decay modes. Fortunately the related 
experimental errors are quite large at present; it is anticipated that further measurements 
with higher precision would observe smaller branching ratios. We also gave the confidence 
levels for some selected decay channels: B^ — > tt'^tt^, K^K~ and B^ — > a;7r+(ir+). 

ACKNOWLEDGMENTS 

We are grateful to S. Laplace for the help on the CKMFitter package. G. Zhu also thanks 
A.I. Sanda, H.n. Li, and Y.Y. Keum for helpful discussions. This work is supported in 
part by National Natural Science Foundation of China. G. Z. thanks JSPS of Japan for 
financial support. 



References 

[1] M. Beneke, G. Buchalla, M. Neubert, and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 

(1999) . 

[2] Y. Y. Keum, H. N. Li, and A. I. Sanda, Phys. Rev. D 63, 054008 (2001); Phys. Lett. 
B 504, 6 (2001); Y.Y. Keum and H.N. Li, Phys. Rev. D 63, 074006 (2001). 

[3] M. Ciuchini, R. Contino, E. Franco, G. MartineUi, and L. Silvestrini, NucL Phys. 
B512, 3 (1998). 

[4] D.S. Du, H.J. Gong, J.F. Sun, D.S. Yang, and G.H. Zhu, Phys. Rev. D 65, 074001 
(2002). 

[5] D.S. Du, H.J. Gong, J.F. Sun, D.S. Yang, and G.H. Zhu, Phys. Rev. D 65, 094025 
(2002); 66, 079904(E) (2002). 

[6] See, for example, A. Hocker, H. Lacker, S. Laplace, and F.Le. Diberder, Eur. Phys. 
J. C 21, 225 (2001); M. Ciuchini et al, J. High Energy Phys. 07, 013 (2001). 

[7] M. Beneke, G. Buchalla, M. Neubert, and C.T. Sachrajda, Nucl. Phys. B606, 245 
(2001). 

[8] BaBar Collaboration, B. Aubert et al, Phys. Rev. Lett. 89, 201802 (2002); Belle 
CoUaboration, K. Abe et al, Phys. Rev. D 66, 071102 (2002). 

[9] For a review, see G. Buchalla, A.J. Buras, and M.E. Lautenbacher, Rev. Mod. Phys. 
68, 1125 (1996). 

[10] J.D. Bjorken, Nucl. Phys. B (Proc. Suppl.) 11, 325 (1989). 

[11] M. Beneke, G. Buchalla, M. Neubert, and C.T. Sachrajda, Nucl. Phys. B591, 313 

(2000) . 

[12] C.W. Bauer, S. Fleming, D. Pirjol, and I.W. Stewart, Phys. Rev. D 63, 114020 

(2001) . 

[13] D.S. Du, D.S. Yang, and G.H. Zhu, Phys. Rev. D 64, 014036 (2001); Phys. Lett. B 
509, 263 (2001). 

[14] M. Neubert and B.D. Pecjak, J. High Energy Phys. 02, 028 (2002). 
[15] M. Beneke, Nucl. Phys. B (Proc. Suppl.) Ill, 62 (2002). 

[16] A. Hocker, H. Lacker, S. Laplace, and F.Le. Diberder, AlP Conf. Proc. 618, 27 

(2002) . 

[17] A. Ah, G. Kramer, and CD. Lii, Phys. Rev. D 58, 094009 (1998). 

[18] M. Neubert and B. Stech, Adv. Ser. Dir. High Energy Phys. 15, 294 (1998). 



[19] T. Feldmann and P. Kroll, Phys. Scripta T 99, 13 (2002). 



[20] P. Ball, J. High Energy Phys. 09, 005 (1998); P. Ball and V.M. Braun, Phys. Rev. 
D 58, 094016 (1998). 

[21] H. Leutwyler, Nucl. Phys. B (Proc. Suppl.) 64, 223 (1998). 

[22] T. Feldmann, P. Kroll, and B. Stech, Phys. Rev. D 58, 114006 (1998); Phys. Lett. B 
449, 339 (1999). 

[23] P. Ball, V.M. Braun, Y. Koike, and K. Tanaka, Nucl. Phys. B529, 323 (1998). 
[24] M. Beneke and T. Feldmann, Nucl. Phys. B592, 3 (2001). 

[25] CLEO Collaboration, D. Cronin-Hennessy et al, Phys. Rev. Lett. 85, 515 (2000); 
CLEO Collaboration, S.J. Richichi et al, Phys. Rev. Lett. 85, 520 (2000); CLEO 
Collaboration, CP. Jessop et a/., Phys. Rev. Lett. 85, 2881 (2000); CLEO Collabo- 
ration, R.A. Briere et al, Phys. Rev. lett. 86, 3718 (2001); CLEO Collaboration, E. 
Eckhart et al, ibid, (to be published), |hep-ex /0206024[ 

[26] BaBar Collaboration, B. Aubert et al, Phys. Rev. Lett, (to be pub- 
lished), pIep^x/0207065; |h^ex/0207055: |hep-ex/0107037| |hep-ex/0107058: 



|hep-ex/ 0109007: Phys. Rev. Lett. 87, 221802 (2001); BaBar Collaboration, A. Bevan, 
talk presented at ICHEP2002, Amsterdam, Holland, 2002; BaBar Collaboration, P. 
Bloom, talk at Meeting of the Division of Particles and Fields, American Physical 
Society, Williamsburg, Virginia, 2002. 

[27] Belle Collaboration, B.C.K. Casey et al, Phys. Rev. D (to be published), 
[hep-ex/0207090; Belle Collaboration, R.S. Lu et al, Phys. Rev. Lett. 89, 191801 
(2002); Belle Collaboration, A. Gordon et al, Phys. Lett. B 542, 183 (2002); Belle 
Collaboration, P. Chang, talk presented at ICHEP2002, Amsterdam, Holland, 2002; 
Belle Collaboration, K.F. Chen, talk presented at ICHEP2002, Amsterdam, Hol- 
land, 2002; Belle Collaboration, H.C. Huang, talk at 37th Rencontres de Moriond on 
Electroweak Interactions and Unified Theories, Les Arcs, France, 2002. 

[28] M. Beneke, talk given at Flavor Physics and CP Violation(FPCP), Philadelphia, 
Pennsylvania, 2002, hep-ph/0207228| 

[29] M. Gronau and J.L. Rosner, Phys. Rev. D 65, 013004 (2002); 65, 079901(E) (2002); 
65, 093012 (2002); 66, 053003 (2002). 

[30] Y.Y. Keum, talk presented at Workshop on CKM Unitarity Triangle (CERN 2002- 



2003), Geneva, Switzerland, 2002, |hi^ph /0209002 Hh^p h/0209208; Y.Y. Keum and 



A.I. Sanda, talk given at 3rd Workshop on Higher Luminosity B Factory, Shonan 
Village, Kanagawa, Japan, 2002, hep-ph/020 9014| 

[31] D.S. Du, C.S. Kim, and Y.D. Yang, Phys. Lett. B 426, 133 (1998); M.R. Ahmady, 
E. Kou, and A. Sugamoto, Phys. Rev. D 58, 014015 (1998); A.L. Kagan and A. A. 
Petrov, hep -ph/970 7354; A.L. Kagan, hep-ph/9806266, 



Y.Y. Keum, |hi^h/0210127[ C.H. Chen, Y.Y. Keum, and H.n. Li, Phys. Rev. D 64, 
112002 (2001). 



